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Foreword 


A  number  of  tables  giving  the  properties  of  air  in 
l-diraensional  flow  have  been  published  in  the  past.  When  a  new 
set  of  tables  such  as  the  present  one  is  issued,  the  question  is 
therefore:  "What  desirable  property  do  these  tables  have  that 
distinguishes  them  from  previous  tables?"  The  answer  is  that  the 
present  tables  are  more  general,  and  expedite  the  solution  of  many 
problems  for  which  conventional  tables  are  not  suited.  The  use 
of  the  Mach  number  as  a  variable  results  in  tables  that  are 
independent  of  the  stagnation  temperature  of  the  gas.  The 
present  tables  go  a  step  farther.  A  new  force  variable  is  intro¬ 
duced  in  such  a  way  as  to  make  the  tables  independent  of 
stagnation  pressure,  while  retaining  the  usual  independence  of 
stagnation  temperature. 

The  project  reference  for  this  work  is  1T061102B11A-07 . 
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ABSTRACT 


New  tables  for  the  compressible  flow  of  air  have  been  prepared, 
based  on  1-dimensional  flow  with  k  -  1.4,  A  new  dimensionless  force 
variable  4  =  kPA/c^m  is  introduced  and  it  is  shown  that  this  variable 
is  independent  of  the  stagnation  pressure  of  the  flow.  This  feature 
makes  the  new  tables  applicable  to  flows  of  different  stagnation  temp¬ 
eratures  and  pressures  in  much  the  same  way  that  conventional  tables 
are  applicable  to  flows  of  different  stagnation  temperatures  only. 

The  argument  of  the  new  tables  is  the  force-variable  6',  the  functions 
tabulated  are  X  (M*),  M,  j,  T/TQ,  c/ca',  P/Po>  f/^o,  A/Aa,  In  (A/Aa), 
and  4fLA/D.  The  tables  are  divided  into  white  regions,  in  which  lin¬ 
ear  interpolation  introduces  no  significant  error;  and  gray  regions, 
in  which  linear  interpolation  can  introduce  error. 


NEW  TABLES  FOR  AIR  FLOW 


1.  Introduction 


New  tables  containing  the  information  needed  for  solving  the 
usual  problems  of  1-dimensional,  compressible  flow  of  air  have  been 
constructed.  Several  advantages  are  claimed  for  them.  They  are 
more  compact  than  previous  tables  yet  they  permit  the  solution  of 
most  of  the  problems  that  can  be  solved  with  more  extensive  tables, 
including  problems  of  general  flow,  isentropic  flow,  Fanno-Iine 
flow,  Rayleigh-line  flow,  the  determination  of  friction  coefficients, 
and  the  calculation  of  properties  across  normal  shocks.-  Many  fluid- 
flow  problems  can  be  solved  more  directly  and  easily  with  the 
present  tables  than  with  previous  tables,  as  will  be  demonstrated 
later  by  examples.  Other  problems  can  be  solved  with  equal  ease 
with  the  present  tables  and  with  previous  tables . 

The  most  important  new  feature  of  the  present  tables  is  the 
use  of  a  new  force -variable'  <( .  The  variable  <f>  has  been  chosen 
because  it  is  easily  computed  from  the  commonly  measured  quantities 
and  because  its  use  greatly  increases  the  generality  of  the  tables 
without  increasing  their  complexity.  Just  as  a  table  employing  M 
rather  than  velocity  can  be  used  for  a  flow  of  any  stagnation 
temperature,  a  table  employing  p  rather  than  pressure  and  area  can 
be  used  for  a  flow  of  any  stagnation  pressure,  as  will  be  demon¬ 
strated  later . 

As  velocity  variables,  both  the  Mach  number  M  and  the 
velocity-number  A  are  given  in  the  present  tables.  The  velocity 
number  X  =  u/c*,  where  u  is  the  stream  velocity  and  c*.  is  the 
velocity  that  the  stream  would  reach  at  a  sonic  throat.  The 
variable A  is  often  represented  by  M  .  In  the  present  report  A  is 
generally  preferred  to  M  because  it  leads  to  simpler  expressions 
for  most  of  the  quantities  to  be  calculated;  M  is  included  for  those 
who  prefer  it  and  for  use  in  those  problems  where  it  is  more  conven¬ 
ient  than  A  . 

The  tables  of  this  report  consist  of  three  small  auxiliary 
tables  and  a  main  table .  The  latter  (Table  4)  contains  all  the 
calculated  functions .  The  functions  are  divided  into  4  groups : 
general,  isoenergetic,  isentropic,  and  frictional.  Each  group  is 
sometimes  referred  to  as  a  table  or  tables  although  all  are 
printed  side  by  side  as  parts  of  Table  4.  Because  </>  is  usually 
the  first  of  the  tabulated  quantities  to  be  found  experimentally, 
it  has  been  chosen  as  the  argument  of  the  table . 


The  first  four  columns  of  Table  4  constitute  the  general  table. 
Most  of  the  new  and  useful  features  of  the  present  tables  are 
associated  with  this  table.  Without  going  into  details  here,  it  will 
be  mentioned  that  this  table  is  used  when  enough  properties  are  known 
to  define  a  state  and  other  properties  of  the  state  are  desired, 
and  that  it  is  also  used  for  other  purposes  by  imposing  various 
conditions  on  it.  Thus,  when  certain  conditions  are  imposed  on  <fi, 
the  general  table  becomes  a  Fanno-line  table. 

Columns  5  and  6  in  Table  4  constitute  the  isoenergetic  table. 

The  reference  states  in  this  table  are  the  stagnation  state 
(subscript  o)  and  the  sonic  state  (subscript*)  .  In  this  table  the 
reference  states  are  isoenergetic  with  the  general  state  under 
consideration.  Stated  otherwise,  the  isoenergetic  tables  apply  in 
situations  where  the  stagnation  temperature  TQ  does  not  change. 

Columns  7  to  10  of  Table  4  constitute  the  isentropic  table.  The 
reference  states  in  this  table  are  the  stagnation  state  and  the  sonic 
state  having  the  same  entropy  as  the  general  state  under  consideration. 
Note  that  in  most  applications  of  the  isentropic  tables  the  flow  is 
both  isentropic  and  isoenergetic.  If  M  is  used  rather  than  the 
present  isentropic  tables  are  entirely  conventional. 

Column  11,  the  final  function  in  Table  4,  constitutes  the 
frictional  table;  it  is  the  conventional  friction  function. 

The  grouping  of  the  tables  in  Table  4,  so  that  the  conditions  oi 
applicability  of  each  group  of  functions  can  be  seen  at  a  glance,  is 
considered  important.  In  previous  tables,  general  relations  have  often 
been  given  in  "isentropic" tables.  Users,  especially  occasional  users, 
have,  often  been  unaware  that  these  relations  were  general.  As  a  result 
many  laborious  and  unnecessary  computations  have  probably  been  made  by 
hand  when  the  tables  could  have  been  used. 

Portions  of  Table  4  are  shaded.  In  the  shaded  portions,  linear 
interpolation  can  result  in  reduced  accuracy.  In  the  unshaded 
portions,  linear  interpolation  will  introduce  no  significant  error. 
Methods  whereby  the  accuracy  of  the  unshaded  portions  can  also  be 
achieved  in  the  shaded  portions  are  presented.  These  methods  Include 
higher-order  interpolation  and  also  the  use  of  the  original  and 
modified  equations. 

Tables  for  the  compressible  flow  of  gases  have  been  published  by 
Keenan  and  Kaye  UL  Emmons  [2  h  the  British  Compressible-Flow  Tables 
Panel  [3],  the  Ames  (NACA)  Research  Staff  {hi,  Bauer  and  Marek[5j,  and 
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many  others.  For  air  the  isentropic  exponent  k  {  =  cp/cv)  is 
usually  taken  as  1.4,  and  this  value  has  been  accepted  in  the 
present  report.  The  present  tables  are  for  air  only.  Conventional 
tables  for  other  values  of  k  are  given  by  Keenan  and  Kaye 
(  k  <>  i.o,  1.1,  1.2,  1.3,  1.4,  and  I.67)  and  by  Bauer  and  Marek. 

The  latter  give  tables  at  intervals  of  0.01  for  k  -  1.01  to 
k  =  1.4l.  The  tables  in  the  present  report  are  intended  for  the 
same  general  purposes  as  those  referred  to  above. 

2 ,  Fluid-Flow  Equations  in  Terms  of  <j>  and  X 

Since  the  variables  used  in  the  present  tables  are  not  the 
conventional  ones,  a  number  of  the  more  important  fluid-flow 
equations  will  be  derived  in  terms  of  the  accepted  variables, 
starting  from  conventional  forms  of  the  equations  of  continuity, 
energy,  and  momentum,  and  using  other  well-known  relations.  Sec¬ 
tions  2,  3,  and  4  of  this  report  give  the  theoretical  basis  of  the 
present  tables  and  other  general  information.  These  sections  may 
be  omitted  by  readers  whose  goal  is  to  learn  how  to  solve  a 
particular  problem.  The  tables,  instructions  for  their  use,  and 
worked  examples  are  given  in  sections  5,  6,  and  7  of  this  report. 

The  new  force  variable  <j  will  first  be  defined;  it  is 
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-  kPA 
c#rfs 


(1) 


where  P  is  the  pressure,  A  is  the  cross-section  area  of  the  flow, 
m  is  the  mass  rate  of  flow,  and  the  quantities  k  and  c*  have  been 
previously  defined.  The  force  is  of  course  PA.  The  normalizing 
factor  k/c^m  makes  <f>  dimensionless  and  gives  it  the  value  1  for 
sonic  flow.  It  will  be  shown  later  that  under  the  usual  conditions 
of  flow,  i>  is  a  function  only  of  A  (or  of  M) .  If  k,  c*,  m,  and  \ 
are  held  constant  and  A  is  varied,  $  does  'not  change.  This  is  true 
because  as  A  increases  P  decreases,  and  vice  versa,  in  such  a  way 
that  the  product  PA  remains  constant.  Stated  otherwise,  is 
independent . of  the  stagnation  pressure  PQ  of  the  flow.  The  very 
useful  fact  that  the  PA-product  is  independent  of  stagnation  pres¬ 
sure  has  not  been  widely  or  clearly  recognized  in  the  past.  Once 
it  is  recognized,  the  advantages  of  using  a  variable  such  as  f 
become  obvious . 


The  constancy  of  the  PA-product  is  made  plausible  by  recall¬ 
ing  some. of  the  concepts  of  kinetic  theory.  According  to  kinetic 
theory,  static  pressure  is  proportional  to  the  rate  at  which 
momentum  is  carried  in  either  direction  across  a  unit  plane  fixed 
relative  to  the  gas,  by  the  thermal  agitation  of  the  molecules. 
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It  turns  out  that  if  we  change  the  cross-section  area  of  a  flow, 
keeping  ifi  and  X  (and  hence  also  T)  constant  we  do  not  change  the 
number  of  molecules  carrying  thermal  momentum  or  the  rate  at  which 
they  carry  it.  In  macroscopic  terms,  when  A  is  doubled  P  is  cut 
in  half  and  the  PA  product  does  not  change. 

Equation  of  Continuity.  The  equation  of  continuity  is 

m  =J>uA  (2) 

where  /»  is  the  density  of  the  flow.  Using  the  ideal  gas  law 
P«J>RT  to  eliminate  f>  ,  the  equation  of  continuity  becomes 

(3) 

ih  u 

We  multiply  this  equation  by  k/c#,  write  T  as  T*(T/T#)  and  u  as 
to  obtain 

kPA  kET#  1  T  (4) 

c*m  c§  A 

The  subscript  star  indicates  as  usual  a  sonic  reference  state. 

In  view  of  the  fact  that 

c  =  qfkR?  (5) 

2 

we  see  that  kRT#/c#  is  equal  to  unity. 

Then,  using  eq  (l),  eq  (4)  reduces  to 

*S  .  i  I  •  (6) 

>.  T„ 


Using  the  relation 


T*r  2 


where  TQ  is  the  stagnation  temperature  of  the  flow,  eq  (6) 
may  also  be  written 


-  tofl  1  T 

This  is  the  desired  form  of  the  equation  of  continuity. 
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Combined  Equations  of  Continuity  and  Energy.  The  energy- 
equation  will  be  accepted  in  the  form 


~T  __  i  \ ^  • 


this  form  is  of  course  valid  only  when  cp,  the  specific  heat  at 
constant  pressure,  is  constant.  When  T/tq  is  eliminated  between 
eqs  (8)  and  (9)  we  obtain 


$ 


k+1  1  „  k-1  . 

2  A  2 


(10) 


which  is  the  desired  form  of  the  combined  equations  of  continuity 
and  energy. 


We  can  now  prove  the  statement  made  earlier  regarding  the 
constancy  of  the  PA-product  of  a  stream.  It  is  clear  from  eq  (10) 
that  is  a  function  of  A  only-,  with  no  restrictions  on  PQ,  T0,  or  m. 
To  examine  the  situation  more  closely  we  will  rewrite  eq  (10), 
restoring  the  value  of  <f>  given  in  eq  (l) 

(11) 

This  equation  shows  that  when  the  quantities  on  the  right-hand  side 
of  the  equation  remain  constant,  P  and  A  can  vary  only  in  such  a  way 
that  their  product  remains  constant. 


The  special  case  in  which  eq  (ll)  is  applied  to  sonic  flow 
( A =  l)  is  of  interest.  In  this  case  we  find 


=  m 

c*  K 


(12) 


The  quantity  P#A#/c^  is  evidently  an  invariant  for  any  1-dimensional 
flow  unless  there  is  addition  or  loss  of  mass,  or  a  process  such  as 
dissociation  that  changes  k.  We  may  change  PQ,  T0,  or  A  in  a  flow  in 
any  physically  possible  way,  changing  P  ,  A*,  and  c#  separately,  but 
the  quantity  P^A^/c*  will  not  change.  In  view  of  eq  (12), 
k/c*m  =  l/p#A#.  and 


^  JSC. 


PA 


(13) 
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This  could  have  been  adopted  as  the  definition  of  instead  of  eq  (l)  * 
The  definition  involving  P#A#  is  perhaps  more  elegant  mathematically, 
but  the  other  is  preferred  for  two  reasons:  first,  eq  (l)  is  more 
convenient  in  practice  because  it  contains  the  experimental  quanti¬ 
ties  from  which  is  usually  calculated;  second,  eq  (l)  avoids  any 
confusion  that  might  arise  concerning  the  meaning  of  P#  and  A*. 

These  quantities  can  have  more  than  one  meaning,  depending  on 
whether  an  isentrope,  a  Fanno  line,  or  a  Rayleigh  line  is  being 
discussed. 

Impulse  Function .  The  impulse  equation,  also  known  as  the 
momentum  equation,  is 

J  =  PA  +  mu.  (1*0 

For  convenience  we  will  refer  to  J  as  the  impulse  of  the  stream. 

Note  that  J  has  the  dimensions  of  force,  whereas  impulse  is  conven¬ 
tionally  the  product  of  force  and  time.  We  multiply  eq  (l4)  by  the  • 
same  factor  previously  used  in 

kJ  _ 
c*™ 

We  define  a  normalized  impulse 

J  * 

so  that  eq  (15)  reduces  to 

0  ' 

Eliminating  <]>  between  eqs  (10) 

(18> 

This  is  the  desired  form  of  the  impulse  equation.  It  is  particularly 
useful  because  it  gives  j  as  a  function  of  A  only,  with  no  reference 
to  P  or  A  or  the  product  PA,  and  with  no  reference  to  T  except  by  way 
of  c#,  which  occurs  in  A.  In  the  derivation  of  eq  (18)  the  equa¬ 
tions  of  continuity,  energy,  and  momentum  were  all  used.  Any 
solution  of  eq  (18)  is  able  to  satisfy  all  three  of  these  equations, 
but  if  information  about  pressures  and  areas  is  required,  it  is 
necessary  to  use  some  other  equation  in  addition  to  eq  (l8) .  For 


the  definition  of  <f>,  obtaining 


kPA  kmu 

+ 


c*m 


c#m 


function 


(15) 


kJ^ 

c*m 


(16) 


kA . 

and  (17)  gives 


1 


(17) 
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example  we  may  know  get  j  from  eq  (l8)j  then  go  to  eq  (17)  to 
find  <fr,  finally  from  <f>  we  get  P  or  whatever  ultimate  quantity  is 
needed. 


The  equations  derived  above  Include  the  most  important 
equations  of  1-dimensional  fluid  flow.  There  are  a  great  many  more 
that  could  be  derived,  but  from  this  point  on  the  treatment  will  be 
confined  mainly  to  those  additional  equations  that  are  actually 
needed  in  connection  with  the  present  tables. 

Equations  (10)  and  (l8)  give  f  and  j  respectively  as  functions 
of  X.  The  inverse  equations  are  often  needed.  Solving  eqs  (10) 
and  (l8)  for  \  we  obtain 


Only  the  positive  sign  of  the  radical  is  retained  in  eq  (19) 
because  A  is  always  positive .  Both  signs  are  retained  in  eq 
(pO);  the +  sign  corresponds  to  supersonic  flow  and  the  -  sign  to 
subsonic  flow.  The  two  conjugate  values  corresponding  to  the  same 
value  of  j  satisfy  the  relation  X q  A2  *  and  correspond  to  the 
velocity  numbers  on  the  two  sides  of  a  normal  shock.  The  minimum 
value  of  j  is  k+1;  lower  values  of  j  are  forbidden. 

Relations  Between  c,  c^,  M,  and  X.  From  eq  (5)  we  obtain 


in  view  of  eq  (6)  there  is  an  alternative  expression 

£  •  (22) 

c# 

The  relation  between  the  Mach  number  M  and  the  velocity 
number  A  can  "be  conveniently  obtained  in  terms  of  c/c^.  Since 
M  =  u/c  and  X  =  u/c# 


From  eqs  (22)  and  (23) a  very  simple  relation  can  be  found; 

n2  =  ^  *  (2^) 

One  of  the  desirable  properties  of  a  mathematical  analysis  is  that 
it  yield  simple,  compact  equations.  Looking  at,  for  example,  eqs 
(6),  10),  (18),  and  (22)  from  this  point  of  view  we  note  that  they 
are  indeed  simple  and  compact.  The  choice  of  <f>  and  \  as  variables 
has  led  to  simple,  convenient  equations  based  on  simple  concepts. 

It  is  possible  that  some  other  choice  of  variables  would  be  still 
better  but  this  does  not  seem  very  likely. 

Isentropic  Relations .  The  isentroplc  relation  between  pres¬ 
sure  and  temperature  is 


the  values  of  T/T  being  obtained  from  eq  (9). 

Similar  relations  based  on  isentroplc  temperature  ratios  hold  for 
other  quantities .  Thus 


?.i) 


w 


However,  when  P/PQ  has  been  found  from  eq  (25)  it  is  more  convenient 

to  obtain  A/A  from  eq  (13)  than  from  eq  (27).  Expressing  P  in  terms 
of  P/P 


Also,  when  A/A^  has  been  found  it  is  easier  to  get  f  jf  from  the 
equation  of  continuity  than  from  eq  (26).  Writing  eq  ^2)  for  a 
general  state  and  for  the  corresponding  isoenergetic  sonic  state  we 
find 

—  =  t±iIo  (29) 

To  A(a/a#)  . 


8 


Stagnation  Pressure .  The  change  in  stagnation  pressure  of 
a  stream  is  a  useful  quantity  closely  related  to  the  entropy  change. 
Changes  are  usually  expressed  as  ratios.  When  values  of  P  are  known 
for  two  states  of  a  flow  the  most  straightforward  way  to  get  the 
stagnation -pressure  ratio  is  to  calculate  j,  get  P/Fq  from  j,  and 
then  use  the  equation 


P 

oa.  _ 


(p/Po>i ' 

(P/Po>2  • 


(30) 


However,  there  is  another  way  of  getting  P  /p  that  is  sometimes 
more  convenient.  oz  oi 

The  alternative  to  eq  (30)  involves  A/A  instead  of  p/p  > 

’*■  o 

and  the  greater  convenience  of  the  alternative  comes  from  the  fact 
that  in  most  problems  areas  remain  constant,  whereas  pressures 
change.  Applying  eq  (12)  to  any  two  states  of  an  isoenergetic  flow  • 
(c  z  constant),  we  see  that 

Tv 

(p#A#)^  _  (p*a*)£,  (31) 

We  solve  this  equation  for  P^2/P#1  and  use  the  fact  that 
p*2/p*i  :  Po2/Pol,  obtaining 


(32) 

poi  A** 


or 


Soar  hk  (A/A*)a  (33) 

poi  A&  (A/A*-)i  . 

If  A1/Ai  remains  constant  during  a  series  of  calculations  and  P2 
does  not,  eq  (33)  is  preferable  to  eq  (30).  Note  that  for  constant 
area  flow  A^/Aj.  drops  out  of  the  equation  . 


Entropy  and  Efficiency.  The  entropy  s  of  a  gas  relative  to 
the  entropy  sr  in  some  reference  state  is 


s-sr 


c 

P 


P 

-R  In  | 

r 


(3M 
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This  equation  is  valid  both  for  flowing  gases  and  for  gases 
at  rest.  If  we  put  s  equal  to  a  constant,  sQ,  eq  (3*0  becomes  the 
equation  of  an  isentrope.  The  point  of  greatest  interest  on  an 
isentrope  is  the  stagnation  state;  its  temperature  is  T0  and  its 
pressure  P0.  Thus 


Sq~s^s  Cp  In  _&  -  R  3n  — o  #  (35) 

-**  r  r  ji 

If  we  use  this  equation  to  calculate  the  change  in  entropy  between 
two  isoenergetic  stagnation  states  (T01<a  T0^)  we  find 

S°2  ~  3q1  =  -  In  Pas.  (36) 

R  Pol  ‘ 

This  quantity  can  be  taken  directly  from  Table  3. 

Friction  Function.  The  friction  factor  f  and  the  friction 
f un c t i orTT+flT/c  are  useful  in  problems  of  1-dimensional  flow  if  f 
is  constant  or  does  not  change  rapidly.  Actual  flows  are  never 
completely  1-dimensional  and  f  is  not  usually  constant,  but  the 
treatment  based  on  these  assumptions  forms  a  convenient  starting 
point  for  more  accurate  treatments  and  is  often  a  good  approxima¬ 
tion  itself.  The  friction  factor  is 


f  5  r 
t  f ui 


(37) 


where  Tis  the  shearing  stress  at  the  wall.  The  force  of  friction 
exerted  on  a  stream  in  a  tube  of  diameter  D  and  length  dl  is 


dF  -  -7 TD  dlT ;  (38) 


using  eq  (37) 


dF  *  -J-1TD  f  f  u^dl  .  (39) 

In  a  constant-area  channel  the  change  in  impulse  dJ  is  equal  to  the 
force  of  friction  dF.  But  from  eqs  (l6)  and  (l8) 


dJ  =  c#m  gL_  -  j^aA  (**0) 

for  an  isoenergetic  flow  (c  =  constant) . 


10 


We  equate  dF  and  dJ  and  eliminate  m  with  eq  (2),  obtaining 


« :  lr-  (s  -  7?)  dA- 


D 


m 


We  will  adopt  the  usual  artifice  of  integrating  this  equation 
along  a  Fanno-line  path,  from  a  general  state  to  the  correspond¬ 
ing  sonic  state,  assuming  a  hypothetical  channel  length  just 
sufficient  to  bring  the  flow  to  the  sonic  state.  Assuming  f 
constant,  the  integral  is 


C\ 
3— - 

k+1 

' 

ln  A  4 

L  u  J 

l  =  2* 

r\  ^ 

A  ’ 


(42) 


The  distance  from  the  general  state  at  £  to  the  sonic  state  at  £ 
will  be  represented  by  L#.  Evaluating  eq  (1*2)  and  multiplying  by 
-2 


-  k.t.1- .  ["2  In  i 
d  “  a  L  A 


(43) 


which  is  the  friction  equation. 


Earlier  in  this  report  it  was  stated  that  the  use  of  A 
usually  led  to  simpler  expressions  than  the  use  of  M.  In  support 
of  this  statement  eq  (43)  will  be  given  in  terms  of  M;  it  is 


kfL^ 

D 


1-M2  k-fl 
kM*  2k  ln 


Cfe+I)  . 

1(1  +'*£rt4)  ’ 


(MO 


eq  (43)  is  considerably  simpler  than  eq  (44). 


Modified  Forms  of  Some  Equations ■  In  some  of  the  gray 
regions  of  Table  4,  especially  those  where  either  ^  or  the  function 
approaches  00,  interpolation  becomes  inconvenient  and  even  impossible 
A  gray- are  a  problem  that  occurs  quite  often  is  to  find  A  from  ?S, 
when  j>  is  large  and  A  is  small.  For  this  situation  it  is  conven¬ 
ient  to  write  eq  (10)  in  the  form 


A 


(45) 
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Thin  equation  must  be  solved  by  successive  approximations  since  the 
unknown,  A  ,  appears  on  both  sides.  However,  when  /  is  large  and  A 
is  small  it  converges  very  rapidly;  it  is  then  easier  to  solve  than 
eq  (l9)*  For  the  first  approximation  the  value  of  A  on  "the  right- 
hand  side  may  be  taken  as  zero,  or  the  nearest  tabular  entry  in 
Table  4  may  be  used. 

Another  way  of  expediting  the  calculation  .of  accurate  values 
In  the  gray  areas  of  Table  4  is  to  use  the  early  terms  of  series 
expansions  of  the  exact  expressions.  As  an  example,  consider  the 
finding  of  P/PQ  when  p  is  large  and  A  is  small.  The  exact  express¬ 
ion  for  P/P  is  eq  (25),  which  is  somewhat  inconvenient  to  evaluate. 
We  eliminate  T/T  from  eq  (25)  with  eq  (9)  and  expand  the  result  by 
the  binomial  theorem  to  obtain 


r  _  ,  k  a  l  k  l  * 
pQ  1  "  Epr  *  +  2  Tc+r~  *  ~  +  •  •  * 


(46) 


This  equation  is  normally  used  only  in  the  region  where  the  terms  in 
X^and  above  can  be  neglected. 

The  basic  equations  in  terms  of  /  and  A  needed  in  the  calcula¬ 
tion  and  the  use  of  the  tables  have  now  been  derived.  A  few  additional 
equations  are  given  without  derivation  in  Section  5  of  this  report. 
Among  them  are  some  series  expansions  of  other  functions  similar  to 
eq  (46) . 


3.  Calculation  of  the  Tables 

The  present  section  describes  the  methods  used  for  calculating 
each -function  given  in  the  present  report.  The  results  are  presented 
in  Section  6  as  figures  (pages  26-28)  and  as  tables  (pages  29-4o) . 
Figures  1  and  2  show  the  various  quantities  as  functions  of  <f>. 

Figure  3  gives  the  same  information  in  the  more  familiar  form  in 
which  M  rather  than  is  the  abscissa. 

Table  1  contains  the  values  of  R  used  in  the  present 

report  and  also  some  often-used  conversion  factors.  Table  2  is  for 
use  in  finding  from  Tq  and  vice  versa.  Writing  eq  (5)  for  sonic 
flow  and  eliminating  T^  with  eq  (7) 
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Table  2  gives  values  of  the  constant  C1  a  -/ 2kR/ (k+l) .  The  need  for 
several  entries  in  the  table  comes  from  the  fact  that  conversion 
factors  have  been  included  in  ,  to  give  c^_  in  any  one  of  three 
sets  of  units,  when  T  is  given  either  in  °R  or  °K. 

Table  3  is  for  use  with  eq  (l) .  For  actual  computations,  it 
is  convenient  to  write  this  equation  in  the  form 


<f 


(49) 


where  C2  is  the  conversion  factor  required  to  reconcile  the  units 
of  PA  with  those  of  c^m,  making  jzS  dimensionless.  Table  3  was  cal¬ 
culated  with  the  aid  of  the  conversion-factor  handbook  of  Zimmerman 
and  Lavine  [6].  Later  the  table  was  verified  by  use  of  the  table 
of  fundamental  constants  and  conversion  factors  given  in  reference  L7J. 
No  discrepancies  were  found. 


In  calculating  Table  4,  the  force -variable  jz 5  was  taken  as  the 
argument  5  A  was  calculated  from  (  using  eq  ( 19)  S  ,1  was  then  found 
from  A  using  eq  (18).  Values  of  T/Tq  were  calculated  from  eq  (9), 
after  which  c/c#  was  found  from 


which  is  a  modified  form  of  eq  (2l) .  Next  M  was  calculated  from 
eq  (23). 

The  isentropic  functions  of  Table  4  were  calculated  from  th= pre¬ 
viously-found  values  of  T/Tq.  First  P/PQ  was  found  from  eq  (25) 
next  A/A#  was  found  from  eq  (28);  finally  f /jo  was  found  from  eq  (29). 
The  friction  function  was  calculated  from  eq  (43) . 

The  length  of  Table  4  represents  a  compromise  between  con¬ 
venience  and  compactness.  The  entries  are  spaced  closely  enough  so 
that  there  are  large  regions  in  which  linear  interpolation  intro¬ 
duces  no  significant  error.  As  mentioned  earlier,  the  white  areas 
of  the  table  permit  linear  interpolation  without  loss  of  accuracy, 
whereas  in  the  areas  with  a  gray  shading  some  loss  of  accuracy  can 
result.  Where  full  accuracy  is  required  in  the  gray  regions,  vari¬ 
ous  alternatives  are  available .  In  some  cases  the  original 
equations  from  which  the  tables  were  calculated  are  convenient  to 
use .  In  others  a  modified  form  of  an  equation,  or  a  series, 
expansion,  may  be  useful.  But  in  the  majority  of  cases  it  will  be 
preferable  to  resort  to  higher-order  interpolation.  Unfortunately 
there  are  regions  where  interpolation  of  any  order  becomes 
impractical.  These  are  the  regions,  generally  at  the  ends  of  the 
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table,  where  either  the  argument  or  the  function  approaches 
infinity.  In  these  regions  some  form  of  the  defining  equations 
should  be  used. 

For  convenience,  the  formulas  for  linear  and  4-point  (cubic) 
interpolation  without  the  use  of  differences  will  be  given.  These 
are  the  formulas  that  are  most  convenient  when  a  desk  calculator  is 
used.  For  linear  interpolation  between  the  two  tabulated  values  f 
and  f,  (do  not  confuse  these  f's  with  the  friction  factor)  the 
formula  is 

f  B  (i-p)fD+pf,  (51) 

where  p  is  the  fraction  of  the  tabular  interval  at  which  f  is 
desired.  Note  that  the  sum  of  the  multipliers  1-p  and  p  is  unity. 
Cumulative  multiplication  should  be  used,  and  if  pf,  is  calculated 
first  the  quantity  1-p  may  be  obtained  as  an  auxiliary  quantity  (the 
counting  dials  are  run  backward) . 

When  higher-order  interpolation  is  used,.  4-point  is  often  used 
rather  than  3-point  (quadratic),  because  4-point  interpolation  is  more 
accurate  and  little,  if  any,  more  difficult.  The  recommended  equation 
is 


f  «.  A_(  f_(  +A0f0  +  A|fj  +Afcfa  .  (52) 

Except  at  the  ends  of  the  table  where  f_,cr  f^may  not  exist,  the  four 
tabular  values  are  chosen  so  that  f  lies  between  fQ  and  f,  .  The  A's 
are  interpolation  coefficients  that  are  functions  of  p,  the  fraction 
of  the  tabular  interval.  Reference  [8]contains  values  of  the  A's 
tabulated  at  intervals  of  0.0001  in  p. 

As  a  check  against  computational  errors,  the  first  and  second 
differences  of  all  functions  in  Table  4  were  machine -computed  and 
printed  as  part  of  the  basic  computer  program.  The  differences, 
especially  were  examined  carefully  for  errors.  Eighteen 
significant  digits  were  carried  in  all  of  the  machine  computations. 
Before  printing  out,  the  functions  were  all  rounded  to  4  decimal  places 
and  the  differences  generally  to  5  decimal  places.  All  rounding  was 
to  the  nearer  figure.  As  an  additional  check  for  errors,  a  number  of 
spot  comparisons  were  made  between  the  present  tables  and  those  df 
Keenan  and  Kayelll  and  the  British  Compressible  Flow  Tables  Panel  [3] • 
The  functions  T/Tq,  P/P  ,  f/f.,  and  A/A  ,  which  appear,  in  all  3  sets 
of  tables,  and  the  functions  \  and  ^)fL  /D,  which  appear  in  reference  Li] 
and  in  the  present  tables,  were  compared  at  M  =  0,  0.1,  0.2...  2.0,  2.5, 
3.0,  4,  5,  10,  and  00. 
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Additional  comparisons  between  reference  J and  the  present  tables 
(at  the  values  of  M  given  above)  were  made  by  using  the  relation 
j  a  (k+l)  (f/f#),  and  by  using  the  fact  that,  along  a  Fanno  line, 

P/P*  is  numerically  equal  to  <f> .  In  reference  [13,  <j  was  obtained 
both  from  the  Fanno-line  table  (Table  42)  and  from  the  isentropic 
table  (Table  30 ) .  For  the  latter,  the  relation  i  -  (P0/P*)  (PA/PQA#) 
was  used . 

The  discrepancies  s hown  up  by  the  spot  comparisons  were 
essentially  those  to  be  expected  from  rounding-off  and  interpola¬ 
tion  errors.  In  the  case  of  a  few  functions  there  were  slight 
indications  of  consistent  discrepancies,  but  these  were  on  the 
border  line  of  detectability.  The  accuracy  of  each  of  the  3  tables 
compared  appears  to  be  adequate  for  all  ordinary  uses. 

The  second  differences  of  the  various  functions  were  used  as 
the  basis  for  dividing  Table  4  into  white  and  gray  areas .  Accord¬ 
ing  to  interpolation  theory  the  maximum  difference  between  the 
results  of  linear  interpolation  and  quadratic  interpolation  reaches 
l/2  unit  in  the  last  decimal  place  of  a  table  when  the  second 
difference  equals  4  in  that  decimal  place.  In  the  white  areas  of 
Table  4  the  second  differences  are  less  than  4  and  in  the  gray  areas 
they  are  normally  4  or  greater. 

4.  Discussion 

The  principal  advantage  claimed  for  the  present  tables  as 
compared  with  previous  tables  is  their  greater  generality.  Whereas 
the  use  of  the  Mach  number  or  velocity  number  makes  a  single  table 
valid  for  flows  of  different  stagnation  temperatures,  the  use  of  the 
force -variable  j)  also  makes  the  same  table  valid  for  flows  of 
different  stagnation  pressures.  The  generality  comes  from  the  way 
in  which  X  (or  M)  and  j  are  defined.  The  tables  are  independent  of 
mass  rate  of  flow  because  of  the  presence  of  ra  in  the  normalizing 
factor  k/c#rfi;  they  are  independent  of  stagnation  temperature 
because  of  the  presence  of  c  in  the  normalizing  factor  and  in  the 
definition  of  X  >  and-  they  are  independent  of  stagnation  pressure 
because  the  factor  k/c^m  is  independent  of  stagnation  pressure. 

Previous  tables  (excluding  shock  tables)  have  usually 
consisted  of  a  single  "isentropic"  table  or  of  an  isentropic  table, 
a  Fanno-line  table,  and  a  Rayleigh-line  table.  The -present  tables 
consist  of  a  general  table,  an  isoenergetic  table,  an  isentropic 
table,  and  a  frictional  table. 
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The  general  table  contains  the  functional  relationships 
that  hold  for  any  1-dimensional  flow.  Since  A  (or  M)  is  uniquely 
determined  by  <3, the  general  table  furnishes  a  convenient  means  of 
finding  X  (or  M)  from  measured  quantities.  Thus  if  we  know  A,  T  , 
and  A  we  can  calculate  c  from  T  and  find  jz$;  then  we  enter  the 
table  to  find  A  *  All  other  tabulated  properties  of  the  flow  are 
also  immediately  available. 


When  the  conditions  of  constant  Tq  (constant  c#)  and 
constant  A  are  imposed  on  <f>y  the  general  table  becomes  equivalent 
to  a  Fanno-line  table.  The  variable  though  defined  in  an 
entirely  independent  way,  is  numerically  equal  to  the  pressure 
ratio  P/P  given  in  a  Fanno-line  table. 

In  Rayleigh-line  flow  the  impulse  J  and  the  area  ofthe  flow 
remain  constant,  but  TQ  (and  c^)  change.  Thus  according  to  eq  (l6) 
the  quantity  c^j  remains  constant  on  a  Rayleigh  line.  By  imposing 
the  conditions  of  constant  A  and  constant  c  j  on  the  general  table 
it  may  be  used  as  a  Rayleigh-line  table .  Tine  conditions  are  not  so 
easy  to  work  with  as  those  involved  in  Fanno-line  flow;  however, 
they  permit  Rayleigh-line  problems  to  be  solved  with  the  general 
table.  In  solving  problems  other  than  Fanno-line  and  Rayleigh- 
line  problems  it  may  be  desirable  to  impose  still  other  conditions 
on  the  general  table .  For  example  we  may  impose  the  condition  of 
constant -temperature  or  of  constant  pressure. 

The  isoenergetic  functions  in  Table  4-  are  separated  from 
the  isentropic  functions  to  emphasize  the  fact  that  an  isoenergetic 
flow  need  not  be  isentropic  and  usually  is  not.  It  is  also  true 
that  an  isentropic  flow  need  not  be  Isoenergetic,  but  it  usually  is. 
The  principal  exception  occurs  when  mechanical  work  is  done  by  or 
on  the  gas.  Thus  in  most  of  the  problems  to  which  the  isentropic 
functions  are  applied  the  flow  is  both  isoenergetic  and  isentropic 
and  both  TQ  and  PQ  are  constant. 

The  quantity  PA/PQA#  tabulated  by  various  authors  is  closely 
related  to  <j  of  the  present  report;  we  see  from  eq  (13)  that 


PA  =  J*  A 
poA*  Po 


(53) 


The  quantity  P  appearing  in  this  equation  must  be  identified  as 
the  local  stagnation  pressure.  If  Po  is  taken  to  be  the  stagna¬ 
tion  pressure  in  some  initial  upstream  state,  and  the  flow  is  not 
isentropic,  eq  (53)  will  not  hold. 
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The  quantity  PA/Pqa#.  has  most  of  the  desirable  properties  of 
f  since  it  is  equal  to  a  constant  (P*/P  )  times  j> .  The  relations 
between  PA/PqA#,  M,  and  F/F*  given  lay  Keenan  and  Kaye  [.l]  and 
by  others  are  general  when  Po  is  interpreted  as  the  local  stagna¬ 
tion  pressure,  although  they  are  presented  as  parts  of  isentropic 
tables*  For  example  the  relation  between  M  and  PA/P  A#  tabulated 
by  Keenan  and  Kaye  may  be  used  to  calculate  M  from  the  pressure 
and  the  flow  per  unit  area,  since  eq  (53)  may  be  written  in  the  form 


PA  . 


kPA 

c*m 


(54) 


After  PA/P  A_x  has  been  calculated  from  this  equation  the  table  is 
entered  and  M  found  directly.  This  is' more  convenient  than  using 
the  mass  flow  parameter  as  described  by  Shapiro  (reference  1 9]  >p82) . 
The  use  of  the  general  section  of  the  present  tables  is  somewhat 
more  convenient  than  either  of  these  methods. 

In  this  report  ^  has  been  used  in  preference  to  PA/PqA^ 
because  of  its  greater  convenience  and  simplicity,  and  because  it 
has  the  value  1  for  sonic  flow  rather  than  the  value  P^/P  .  Note 
however  that,  when  is  less  than  1,  X  and  M  are  greater  than  1,  and 
vice  versa. 


Since  a  high  degree  of  generality  is  claimed  for  the  present 
tables,  mention  should  be  made  of  the  generalized  compressible-flow 
function  of  Benedict  and  Steltz JloJ .  This  function  is 


(55) 


the  authors  show  that  it  satisfies  the  equation 
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(56) 


With  eq  (56)  the  generality  of  the  present  treatment  is  achieved  in  a 
somewhat  roundabout  way.  It  can  easily  be  shown  that  P  is  not  really 
a  new  function  at  all,  that  in  fact 


(57) 


where  A  /A  is  the  isentropic  area-relation. 

The  rather  lengthy  treatment  of  critical  states  and  of 
entropies  given  by  Benedict  and  Steltz  is  unnecessary,  as  was  pointed 
out  by  A.  J.  W.  Smith  in  his  discussion  of  the  paper  (reference  £l0], 
p  68) .  Smith  derived  a  function  X,  previously  used  by  others  in 
England,  and  showed  that  p  was  equal  to  a  constant  times  X.  It  is  not 
clear  whether  any  of  these  workers  were  aware  that  their  functions 
were  simply  variations  of  A  /A.  Probably  Benedict  and  Steltz 
considered  T  as  a  function  £hat  in  some  situations  was  equal  to  kjk 
but  in  other  situations  had  a  different  meaning.  At  any  rate,  they 
did  not  emphasize  the  relationship,  nor  point  out  to  the  reader  that 
the  reciprocal  of  the  function  P  had  already  been  given  in  existing 
tables . 


It  will  be  shown  that  eq  (56)  actually  introduces  the  independ¬ 
ence  of  stagnation  pressures  that  is  achieved  in  the  present  treatment 
by  the  use  of  <f>.  We  first  eliminate  and  rj,  from  eq  (56)  by  using 

eq  (57),  obtaining 


(58) 


The  first  term  in  this  equation  may  be  written  as  c^/c  ,  and  the 
second  as  P^/P*,  >  giving  '  ’ 


or 


C  *1 

P*a 

A*a 

c*a 

p*. 

A*l 
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(59) 


Since  this  relation  holds  between  any  two  points  in  a  flow,  the 
quantity  P^A^/c^  must  be  a  constant  for  the  flow.  We  do  not  get  the 
value  of  this  constant  from  the  equation  of  Benedict  and  Steltz,  but 
from  eq  (12)  we  see  that  it  is  m/k.  Since  P/P#  and  A/A^  are  known  to 
be  functions  of  X  (or  M)  only,  PA/c  must  likewise  be  a  function  of  A 
only,  and  the  generality  of  eq  (56)  has  been  demonstrated. 

Approximate  solutions  to  difficult  problems  can  often  be 
obtained  by  assuming  the  flow  to  be  isentropic.  In  other  cases  the 
solution  to  a  problem  can  be  obtained  for  the  isentropic  case,  but 
the  result  actually  is  not  limited  to  isentropic  flow.  Benedict  and 
Steltz  derived  eq  (5 6),  which  contains  a  number  of  isentropic  relations, 
but,  as  they  point  out,  their  result  is  not  limited  to  isentropic  flow. 
Other  authors  have  also  developed  "Isentropic"  functions  that  are  not  . 
limited  to  isentropic  flow.  When  practical  it  is  desirable  to  derive 
general  relations  first,  and  limit  them  to  special  cases  only  when 
necessary.  This  has  been  done  in  the  present  treatment. 

5.  Instructions  and  Working  Equations 

This  section  of  the  present  report  contains  the  rules  to  be  followed 
and  the  equations  to  be  used  in  solving  problems .  Much  of  the  infor¬ 
mation  is  taken  from  earlier  sections  of  the  report  and  generally  no 
attempt  is  made  here  to  give  full  explanations.  The  working  equations 
are  all  grouped  together  at  the  end  of  this  section.  A  number  of 
equations  not  given  in  the  main  text  of  this  report  have  been  included 
for  completeness.  The  equations  in  this  section  are  all  designated  by 
letters:  (a),  (b),  (c),  etc.,  and  all  references  to  equations  in  this 
section  of  the  report  are  made  in  terms  of  these  letters. 

The  instructions  are  given  in  the  form  most  suitable  for  finding 
other  quantities  when  m,  A,  T0,  and  P  are  known.  In  many  situations 
some  of  these  quantities  will  be  unknown  and  other  quantities  will  be 
known.  The  user  of  the  tables  should  in  such  instances  have  no  great 
difficulty  in  inverting  or  otherwise  revising  the  instructions  given. 

To  Calculate  The  first  operation  to  be  performed  in  using 
the  tables  is  usually  to  find  the  force -variable  ^j.as.the  first  step 
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in  finding  the  sonic  velocity  c  must  be  found.  Equation  (c)  is 
used  for  this  purpose,  with  a  value  of  the  constant  Cj  taken  from 
Table  2. 

When  c*  has  been  found,  the  force-variable  <f>  is  calculated 
from  eg.  (e),  using  the  appropriate  value  of  taken  from  Table  3* 

Xn  many  applications  it  will  be  found  that  ^  is  egual  to  a  constant 
times  P  and  computations  are  expedited  by  first  finding  the  constant 
and  then  multiplying  by  all  the  various  values  of  P. 

Using  Table  4,  Interpolation .  When  <f>  has  been  found,  the 
corresponding  values  of  X,  M,  j,  and  all  the  other  tabulated 
functions  can  be  obtained  from  Table  4.  Linear  interpolation  is  to 
be  used  in  all  white  areas  of  the  table.  In  these  areas  the  second 
differences  of  the  tabulated  functions  are  less  than  4  in  the  last 
decimal  place  and  the  errors  introduced  by  interpolation  may  be 
expected  to  be  less  than  l/2  in  this  place.  For  linear  interpolation, 
eg  (av)  should  be  used.  In  the  gray  areas  of  Table  4,  second 
differences  are  ^  4  and  linear  interpolation  can  introduce  errors  in 
the  last  decimal  place  or  even  preceding  it.  In  the  gray  areas,  4- 
point  interpolation  is  often  satisfactory;  for  4-point  interpolation 
eg  (aw)  should  be  used.  To  verify  the  accuracy  of  4-point  interpo¬ 
lation,  5-point  interpolation  may  be  used.  If  the  4-point  and  5" 
point  interpolations  agree,  it  may  be  assumed  that  no  error  is  being 
introduced. 

Alternatives  to  Interpolation.  In  many  instances  the 
procedures  given  in  the  table  below  are  preferable  to  4-point 
interpolation,  being  more  convenient,  more  accurate,  or  both. 


Desired 

Recommended  Method 

Function 

for  Gray  Areas 

X 

eg  (Jl) 

J 

eg  (t) • 

T/T 

eg  (v) 

c/°* 

eg  (w) 

P/po 

eg  (ac) 

f/fo 

eg  (ag)  or  (ah) 

V*. 

eg  (am) 
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Of  course  a  function  can  always  be  calculated  to  any  desired 
accuracy  from  the  equations  used  in  the  preparation  of  Table  4, 
Equations  (t),  (v),  and  (am)  were  in  fact  used  in  preparing  Table  4, 
Unfortunately,  many  of  the  other  primary  equations  are  inconvenient 
to 'work  with. 

Other  Operations .  To  find  stagnation-pressure  ratios 
(PQa/P  ),  use  either  eq  (ap)  or  eq  (aq) .  The  latter  is  more 
convenient  when  A(/Aj,  is  constant.  When  In  (Poa/PQ,)  Is  required, 
use  eq  (ar),  getting  the  values  of  In  (A/A  )  from  Table  4. 

To  solve  Fanno-line  or  Rayleigh-line  problems,  use  the  general 
table  subject  to  the  restriction  of  eq  (f)  or  eq  (u),  as  the  case  may 
be . 


To  solve  normal-shock  problems,  first  find  the  desired  proper¬ 
ties  of  the  flow  before  shock  from  Table  4-,  Then  using  eq  (at),  find 
the  value  of  A  after  a  normal  shock.  Enter  Table  4  at  this  value  of 
X  and  find  the  desired  properties  of  the  flow  after  the  shock.  Note 
that  eqs  (at)  and  (au)  are  equivalent  and  either  one  may  be  used. 

Equations .  The  equations  given  below  consist  of  the  more  use¬ 
ful  ones  from  the  main  text  plus  a  number  of  other  useful  equations 
that  were  not  needed  in  the  main  text.  The  latter  group  of 
equations  are  designated  only  by  letters.  Where  both  a  letter  and 
a  number  appear,  the  number  is  the  one  used  in  the  main  text. 


c  =  -/kRT 

(a) 

(5) 

c„=  /2k  R  T 

(b) 

(47) 

w  y  K£L  ° 

V  c,VF 

(c) 

(48) 

i  *  s*. 

(a) 

(1) 

c#ra 

$  -  C  kEA 

(a) 

(49) 

c  ra 

^  =  constant  x  P 

(f) 

(Fanno-line  flow) 
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i 

* 


k+l  1  k-l 
2  h~  2 

k+l  1  T 
2  >  Tq 


A 


FA 


A  * 

A  = 


k+l 

“5“ 


_J—  + 

k+l  “ 


* 

m  =  j°uA 

m  o  EM 
RT 


Ms  U 

c 


■(g)  do) 

00  (8) 

(l)  (13) 

U) 

(k)  (19) 

(A)  (*»5) 

(m)  (20) 

(n)  (2) 

(o)  (3) 

(p) 

(q)  (23) 
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J  = 

c*m 

( r ) 

(16) 

<j  s  ^  +k  X 

(s) 

(17) 

*  ■  ¥(*+-*) 

(t) 

(18) 

c#J =  constant 
(Rayleigh -line  flow) 

(u) 

J-  =  1  -  54  A1 

T  k+1 

0 

(v) 

(9) 

£  .fiT 
c* 

(w) 

(22) 

c  =  pf" 

%  1~* 

(*> 

(21) 

C  J~^Q.  R. 

«.  7M*o 

(y) 

(50) 

•  1.095^ 

(k-1.4) 

(*) 

k 

P  .  /  T  \  k-1 

fo  "  l  ToJ 

(aa) 

(65) 

P  r  1--*  A2+  A  -* 

P„  k-fl  2  k+1 

-*•  _  -+ 
k+1  A 

(ab) 

va 

=  1  -  0.5833  >  '  (ac) 

(k  =  l.h,  error <5  x  10~5  for  0<  \  ^  0.1^2) 
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In  JL  * 


k 

k+-l 


*-1  \4+ 

E£T  A  +  3  k+l 


JL_  Jf  ,  1  2~K 
k+lA  *"  2  (k+l)* 


=  i-0 .4167  X2 

(&a  1.4,  error -^5  x  10  ^  for  0.<)^<0 .176) 


(ad) 


(ae)  (26) 

(af) 


(ag) 


t  «  P/Eo- 

f°  T/T0 

f  f »/  fa 

f0"  A(A/A*T 

L.  it  j.  /.l' 

Aw  fo  A  \  T0 


A_  -  p*  jL 

\  p  p/p 

^  o  '  o 


=  0.5283  _4-  “  !•**) 

P/Fo 

.  gig  in  A  +  j,  -  I  ] 

{  =  -  A(4fL *M 

A (4l/D) 


(ah) 

(ai)  (29) 

(aj) 

(ak) 

(al)  (28) 

(am) 

(an)  (43) 

(ao)  ■ 


24 


(ap)  (30) 


P°»  - 

Poi  P,  (P/P0)g 

P0a.  A,(A/A#)a 

Pot  Aa(A/A#)i 


(aq)  (33) 


■SoC-<?.tS  -l.n  -Z°*  (TQ  :  constant)  <as)  ^6) 

»  poi 

»  „  jr .  (normal  shock)  (a^) 

V  A, 

j  -  (normal  shock)  (au) 


f  s  (l-p)fD+pf,  (av)  (51) 

(linear  Interpolation,  p  is  fraction  of  interval) 


f  =  A.,  fw+Ao^A.  VAzfa.  ,  „  fQl  v 

(4 -point  interpolation,  take  A's  from  reference  L8J.) 


(52) 


’6.  Figures  and  Tables 

The  new  functions  are  given  in  the  present  section  of  this 
report.  They  are  shown  graphically  in  Figs.  1,  2,  and  3,  and  are 
given  in  tabular  form  in  Table  4.  Figures  1  and  2  show  the  various 
functions  with  the  new  variable  $  as  abscissa.  Figure  3  gives  the 
same  information  in  a  more  conventional  form,  with  M  as  abscissa. 

Tables  1,  2,  and  3  are  auxiliary  tables  containing  various  constants 
and  conversion  factors  that  are  often  needed  when  the  main  table 
(Table  4)  is  being  used. 

The  methods  used  in  calculating  the  functions  are  given  in 
Section  3,  and  examples  of  the  use  of  the  tables  are  given  In  Section  7  • 
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Fig.  2.  Four  of  the  tabulated  functions  (T/T q,  c/c¥_)  P/P and- 
f /f*0)  plotted  versus  <{>. 
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©w  V# 

.05  .1  .5  1  5  10 


M  — — -s- 

Fig.  3*  Nine  of  the  tabulated  functions  (X,  j,  T/T  >  c/c  ,  P/P  , 
f/fo'  A/A^,  ln(A/A^),  and  4fL  /d)  and  ^  pl8tted  versus°M. 
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Table  1.  Constants  and  conversion  factors  . 


H  (all  gases)  -  8.31^3  J/mol  K 
mol  wt  (air)  =  28.970  g/mol 
R  (air)  =  0.28700  j/g°K 

o 

s  0.06859^  thermochemical  cal/g  K 

o 

a  0.068594  thermochemical  Btu/lbm  R 
,  o 

=  53 . 3*4-2  ft  lbf/lbm  R 

=  1716.2  ft2/secgoR 

1  Btu/lbm  r  25/020  ft  /  sec 
_L 

1  (Btu/lbm)4 =  158.18  ft/sec 

~  48.213  m/sec 

1  atm  s  760  mm  Hg  (Torr) 

=  14.696  psia 

1  psia=  51*715  Torr 

T  (°H).  t  (°R)  +  459.67 

-  1.8T  (°K) 
o 

T(  K)  s  t  (  C) +  273.15 

t  (°F)  +  40  r  1.8  [t  (°C)+4oJ 


^Except  for  the  molecular  weight  of  air,  which  was  taken  from 
reference  [1],  the  values  in  this  table  were  taken  or  calculated 
from  reference  [ 7 j,  which  uses  the  physical  constants  approved  by  the 
National  Academy  of  Sciences  and  the  National  Bureau  of  Standards. 


Table  2.  Values  of  C,  for  use  in  the  equation  c# 
c#  and  Tq  are  expressed  in  the  indicated  units. 


when 


T 

\  0 

c  \ 

*  \ 

°R 

(=  °F  abs) 

°K 

/Btu\i 

\lbmj 

O.28289 

0.3795^ 

ft 

sec 

I4.4 . 7^7 

60.034 

m 

sec 

13.639 

18.298 

Table  3.  Values  of  C2 required  to  make  ^  dimensionless  in  the 
equation 

Find  the  row  headed  by  the  dimensions  of  PA  and  the  column  headed 
by  the  dimensions  of  c#m  and  use  the  value  of  C^lying  in  the  row 
and  column  thus  found. 


kPA 


PA  \ 

cm  g 

sec  sec 

ft  lbm 
sec  sec 

ft  lbm 
sec  hr 

/  Btui^lbm 
\lhry  sec 

1 

/Btu^lbm 

hr 

dyne 

1 

7.2330 
x  10t5 

0.26039 

4.5727 
x  io“? 

1.6462 

x  10 

lbf 

4.4402 

5 

x  10 

■3EBH 

1.1583 

5 

X  10 

0.20340 

732.26 

atm  cm^ 

1.0132 

x  - 

73.289 

2.6384 

5 

x  10^ 

0.46333 

1668,0 

atm  in3 

■at-kw 

MM 

472.83 

1.7022 
x  106 

2.9892 

10,761 

mm  Hgcm2 

1333.2 

0.096432 

347.16 

6.0965 

x  10 

2.1947 

mm  Hg  in 

8801.4 

0.62214 

2239.7 

3.9332 

-3 

x  10 

14.159 
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Tattle  1*.  Functions  for  the  1-dimensional  flow  of  air  (k=lA). 
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Table  4.  Functions  for  the  I- dimensional  flow  of  air  (k=1.4)J  continued. 


L3 

B 
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its*  to  04  T-l 
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2601  0.3822  1,198?  >  9. IS 09  0,1476 


Table  4.  Functions  for  the  1-dimensional  flow  of  air  (k=1.4),  continued. 
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Table  4.  Functions  for  the  1- dimensional  flow  of  air  (k=1.4),  continued. 
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Table  k.  Functions  for  the  1-dinehsional  flow  of  air  (k='i . k ) ,  continued 
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Table  4.  Functions  for  the  1-dimensional  flow  of  air  (k=l.l;),  continued. 
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Table  k.  Functions  for  the  1-dimensional  flow  of  air  (h=1.4),  continued. 
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Table  4.  Functions  for  the  l-dimensional  flow  of  air  (k=1.4),  continued. 
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Table  4.  Functions  for  the  1- dimensional  flow  of  air  (k=1.4),  continued. 
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Table  4.  Functions  for  the  1-dimensional  flow  of  air  (k=1.4),  continued 
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7 .  Examples  of  the  Use  of  the  Tables 


This  section  contains  9  illustrative  examples  of  the  use  of 
the  present  tables .  After  a  problem  has  been  stated,  the  method  of 
solution  is  described.  Then  the  given  numerical  values  are  listed, 
and  finally  the  computed  values  are  given,  in  the  order  in  which 
they  are  obtained.  The  earlier  examples  are  simpler  than  the  later 
ones,  and  are  explained  in  greater  detail.  For  example,  the  methods 
of  finding  e  and  </>  are  described  only  in  problem  1,  although  these 
quantities  are  used  in  most  of  the  other  problems  also.  The  • 
equations  referred  to  in  the  examples  will  be  found  in  Section  5 
of  this  report. 

Example  1.  Finding  M  _or  \  (also  finding  c*  and  b) .  The  stag¬ 
nation  temperature,  pressure,  flow  rate,  and  cross-section  area  of  a 
stream  of  air  are  known.  Find  the  velocity  number  \  and  the  Mach 
number  M  of  the  flow. 

Solution.  First  find  using  eq  (c)  and  Table  2.  Then  find  i 
using  eq  (e  )  and  Table  3.  Then  enter  Table  4  with  to  find  \  ,  M, 
and  any  other  desired  properties. 


Given: 

T  =  550°R 
o 

P  =6  psia 
m  =0.5  lbm/sec 
A  =11  in2 


Find : 

c, 

i 

i> 

X 


=  0.28289  tTt~ 

*  o 

=  6.63*t  Btu1,/2/lbm1//2 

_  r  kPA 
“  2 

=  0.20340  kPA 


=  5 • 666 
=  0.2103 
=  0,192? 


/ 


hi 


Example  2.  Finding  other  properties  when  M  _is  given.  The  Mach 
number  of  a  flow  of  air  is  known  from  the  angle  of  a  weak  shock  wave 
observed  in  the  flow.  The  stagnation  temperature  and  the  static 
pressure  of  the  flow  are  known.  Find  the  temperature,  stagnation 
pressure,  and  flow  per  unit  area. 


Solution .  Enter  Table  4  with  M  to  find  si.  Use  linear  interpola¬ 
tion  first  but,  since  M  lies  in  a  gray  area,  check  the  result  by  4- 
point  interpolation.  A  small  correction  in  <j>  is  required  to  give  the 
observed  value  of  M  exactly.  Since  the  difference  between  the  "linear" 
and  "A-point"  values  of  i>  is  small,  it  is  safe  to  accept  the  4-point 
value  as  correct.  Using  this  value  of  <(>,  find  T/T  and  P/P  from 
Table  4;  from  them  find  T  and  P  .  Note  that  T/T  and  P/P  iie  in  white 
areas,  so  linear  interpolation  is  adequate.  Nex?  solve  eq  (  e)  for 
m/A,  and  note  that  all  quantities  required  for  finding  m/A  are  known 
except  c.A,  which  can  be  obtained  from  . 


Given: 


M  =  2.6 

T  =  560°R 
o 

P  -2  psia 


Find: 

si  (linear) 

&  (4-point) 

T/T 

P/P° 

T  °  ' 

P 

o 

♦ 

m 

A 

c* 

m 

A 


0.2748 
0.2747 
0.4251 
0.0501 
238. 1°R 
39.9  psia 

c  JsL_ 

c*<i 

1059  ft/sec 

2 

0.3097  lbm/in  sec 


Example  3.  Choked  isentropic  flow,  finding  m.  Air  flows  isen- 
tropically  from  a  tank  where  its  stagnation  temperature  and  pressure 
are  known  into  a  converging  nozzle  of  known  throat  area.  The  flow  at 
the  throat  is  sonic  (choked).  Find  the  general  formula  for  the  mass 
rate  of  flow  in  terms  of  stagnation  quantities  and  throat  area,  and 
apply  it  to  the  given  problem. 
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Solution.  Using  eq  (  e)  and  eliminating  cA  with  eq  (c),  a  general 
expression  is  obtained  that  may  be  solved  for  m.  Writing  this  expres¬ 
sion  for  sonic  flow  =  1) 


« 

in 


^2 


Since  is  known  rather  than  this  may  be  written 


m 


=  C2 


P  A. 
o  * 


? 


This  is  the  general  formula. 


To  apply  the  general  formula  to  a  particular  problem,  the  units 
of  the  quantities  involved  must  be  specified.  Then,  getting  C.,  from 
Table  2,  C  from  Table  3,  and  P^/P  from  Table  4,  the  given  values  are 
substituted  and  m  is  found. 


Given: 

T  =  6G0°R 
o 

Pq  =  20  psia 
Aa  =0.1  in2 

Accept: 

lbf  as  force  unit 
ft/sec  as  velocity  unit 
lbm/hr  as  flow  unit 


Find: 

Cx  =  44.747 

C2  =  1.1583  x  105 

P./P  =  0.5283 
*  o 

m  -  1914  PoA* 


^3 


=  156.3  lbm/hr 


Example  4,  Isentropic  flow,  unchoked.  Air  flows  isentropically 
in  a  converging  nozzle,  but  the  flow  at  the  smallest  cross-section 
(exit)  is  not  choked.  The  stagnation  temperature  and  pressure  of  the 
flow  are  known;  and  also  the  flow  per  unit  area  at  the  exit.  Find 
the  other  properties  of  state  e,  the  exit  state. 

Solution.  First  find  m/AA,  for  the  sonic  reference  state  of  the 
isentropic  flow.  For  this  state  <t>  =  1  and  eq  (  e  )  gives  m/A^  =  C2kPA/ 
cA.  Finding  PA  from  P0  and  c  from  T0  we  calculate  m/A^.  Since  1T1/ A^ 
is  given,  we  can  calculate  A  / A^.  This  permits  us  to  enter  Table  4, 
where  all  the  desired  properties  of  state  e  may  be  obtained.  The 
flow  is  assumed  to  be  subsonic.  Note  that  any  value  of  A/A^  occurs 
twice  in  Table  4,  once  in  the  subsonic  range  and  once  in  the  super¬ 
sonic  range. 

Given; 

P0  =  30  psia 

To  -  1000 “a  2 

m/A  =  0.35  lbm/in  sec 
e 


=  32.174 
=  15.85  psia 
-  44.747 
=  1415  ft/sec 

=  0.5045  lbm/sec 


=  1.441 


=  2.370 

P^PQ  =  0.8687 
=  26.06  psia 

M  =0. 4530 
e 

Example  5.  Fanno-line  flow.  A  flow  of  air  with  friction  at 
constant  area  and  T0  (Fanno-line  flow)  has  pressure  P^  at  station  1 
and  is  sonic  at  the  channel  exit.  The  stagnation  temperature,  the 
mass  rate,  of  flow  and  the  cross-section  area  of  the  channel  are 
known.  Find  the  pressure  at  the  channel  exit  and  the  velocity  of 
flow  at  station  1  and  at  the  exit. 


Find; 


m 

A, 


Solution.  In  Fanno-line  flow,  $5  is  proportional  to  the  pressure. 
We  first  find  the  constant  of  proportionality,  calculating  c.  from  TQ 
and  evaluating  the  parentheses  in  the  equation  ^  =  (G^kA/c^ifi)  P. 

Then,  putting  P  =  P1  we  find  ^  from  which  At  may  be  found.  To  find 
u-,  we  use  eq  (j).  The  pressure  PA  at  the  exit  is  found  by  putting 
&  —  1  in  the  equation  relating  ^  and  P.  The  velocity  at  the  exit  is 
simply  cA. 

Given: 

T0  =  530°R 

m  =  650  lbm/hr 
2 

A  =0.25  in 
P^  =  31  psia 

Find: 

=  44.747 

cA  =  1030  ft/sec 

c2  =  1.1583  X  105 

C^hA/c^m  =  0.0605H 

&  =  0.06051#  (P  in  psia) 

=  1.877 

=  0-600  9 

u^  =  619-1  ft/sec 
P*  =  16.52  psia 


Example  6.  Rayleigh-line  flow.  Air  flows  in  a  constant-area 
channel  without  friction,  but  with  heat  addition.  The  velocity  num¬ 
ber  and  the  stagnation  temperature  are  known  at  station  1.  Down¬ 
stream,  heat  is  added  and  the  flow  chokes  at  the  channel  exit,  sta¬ 
tion  2.  Find  the  stagnation  temperature  when  choking  occurs. 


Solution.  In  a  Rayleigh-line  flow  the  impulse  J  remains  con¬ 
stant.  Since  j  =  kJ/cAm,  the  quantity  cAj  remains  constant  when  J 
is  constant.  Entering  Table  4  with  Al  we  f ind  j 1 . 

ft  1  .■#  J-i 


From  T„.  we 


-  °1 

calculate  c^j..  Then  cAj  for  the  Rayleigh-line  can  be  found.  Since 
the  exit  flow  is  choked,  j  ^  is  known.  Equating  "the  known 

value  c  j  ,  we  can'  find  c  and  from  c  "  we  then  "get  T 

#1  1  •* 


#2 


*a 


02 


t 


Given: 


X1  =  0.6 

T  =  530°R 

=1  (choked  flow) 

Find: 

C  =  44.747 

c.-  =  1030  ft/sec 
*1 

j  =  2.7200 

C*l^l  =  2802  ft^sec 

j2  =  2.400 

cA2  =  1168  ft/sec 

T  _  =  681° R 
o2 

Example  7.  Isothermal  flow  and  choking.  Air  flows  isothermally 
with  friction  in  a  constant— area  channel.  Heat  exchange  between  the 
stream  and  its  surroundings  takes  place  by  some  mechanism  (such  as 
radiation)  that  avoids  the  coupling  that  usually  exists  between  heat 
transfer  and  friction.  The  temperature,  flow  per  unit  area,  and  in¬ 
let  (state  1)  pressure  are  known.  The  downstream  pressure  is  suffi¬ 
ciently  low  so  that  the  flow  at  the  channel  exit  (state  2)  is  choked. 
Find  the  stagnation  temperatures  and  pressures  of  both  states  and  the 
static  pressure  of  state  2. 

Solution.  When  T  is  known  and  is  not,  it  is  easier  to  find  M 
than  either  ^  or  ^  .  First  calculate  the' stream  velocity^at  the  in¬ 
let  (u-j )  using  the  equation  of  continuity  in  the  form  u  =  mRT/PA,  and 
taking  rrom  Table  1  the  value  of  R  that  will  give  u  in  the  desired 
units.  Next  calculate  c  =  /kRT  taking  from  Table  1  the  value  of  R 
that  will  give  c  in  the  same  units  as  those  of  u.  Calculate  s  u^/ c 

enter  Table  4  with  M  .  Find  then  T/To;,  ,  which  gives  T0l;  and 
P  /p  which  gives  F0,.  To  calculate  state  2,  make  use  of  the  fact 
that  isothermal  choking  occurs  when  M  =  1/  J  k.  Entering  Table  4  with 

this  value  of  M^,  find  «L»  T/ ^ oo*  anc*  ^2^°2*  From  T/T02  §efc  T°2* 

From  eq  (e),  P  =  ^c^m/C^kA.  Get  c.„  from  and  C2  from  Table  3 

and  calculate  Finally,  with  P2/P02  calculate  P02. 


1+6 


Given: 


T  =  1200°R  2 

m/A  =  0,25  lbm/in  sec 

=  40  psia 

M2  =  l//k  =  0.8452 

Find: 

R  =  53.342  ft  lbf/lbm°R 
400.1  ft/sec 

R  =  1716.2  ft'^/sec^R 
c  =  1698  ft/sec 
M  =  0.2356 

4.624 

T/T0  ^0.9890 
T01  =  1213°R 
P1/Po1  =  0.9621 
PQ1  =  41.58  psia 
\  -  1.212 
T/T0  2  =  0.8750 
T0  2  =  1371°R 
P2  -  <2s^c2» 
c.„  =  1657  ft/sec 
C2  =  32.174 
?2  =  11.15  psia 

VPo2  -  0.6266 

P02  -  17.80  psia 


47 


Example  8.  Friction.  An  air  flow  is  carried  in  a  circular 
channel  of  constant  cross  section;  it  is  isoenergetic  but  subject 
to  friction.  The  stagnation  temperature  and  flow  per  unit  area  are 
known.  The  pressure  is  measured  at  station  1  and  at  station  2,  the 
latter  being  a  known  number  of  pipe  diameters  downstream  from  the 
former.  Find  the  friction  factor  f  ,  assuming  it  to  be  constant 
between  the  two  stations. 

Solution.  The  friction  equation  (aa)  must  be  evaluated.  First 
calculate  and  {$7;  enter  Table  4  with  these  values  and  find 
(4fLA/D),  and  C4^>Cjfe/D) ^ -  Since  the  friction  function  lies  in  a  gray 
area,  both  linear  and  4-point  interpolation  are  used.  Corresponding 
values  are  . so  near  together  that  the  4-point  results  are  probably 
mathematically  correct.  However,  in  view  of  the  relatively  low  ac¬ 
curacy  of  most  friction  calculations,  the  values  of  4fLA/D  are  rounded 
off  to  3  decimal  places.  The  value  of  f  is  calculated  from  eq  (ao). 


Given: 

T0  =  700°R 

m/A  =  830  lbm/hr  in2 


-  18  psia 
P2  =  16  psia 


A (L/D)  =  60 

Find: 

cA  =  .1184  ft/sec 


=  2.970 
i2  =  2.640 

(4fLA/D)1  =  3.0774  linear 


=  3.0770 
=  3.077 

(4fL*/D)2  =  2.1649 
=  2.1647 
=  2.165 

A(4fL*/D)  =  0.912 
f  =  0.00380 


4-point 
rounded 
linear 
4-po int 
rounded 
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Example  9.  Normal  shock.  A  flow  of  air  undergoes  a  normal 
shock.  Given  the  upstream  values  of  P  and  4,  find  the  pressure 
after  the  shock;  also  find  the  stagnation  pressures  before  and  after 
the  shock. 


Solution,  We  use  the  fact  that  across  a  normal  shock  =  1  /A, 
eq  (at).  Entering  Table  4  with  4,  we  get  A ■  Calculating  \2  1 

from  A ^  we  next  get  ?S-  from  Table  4.  Referring  to  eq  (id)  we  see 
that  all  quantities  in  it  remain  constant  across  a  normal  shock  ex¬ 
cept  4  and  P.  Hence  =  ^2^1  an<^  ^2  C3n  ^oun<^‘  Then,  getting 

and  ^2^02  from  Taale  47  we  calculate  PQ^  and  P02. 

Given: 

P^  =  4  psia 

41  =  0.5 

Find: 

\x  =  1-5 
X2  =  0.6667 

4 2  -  1.6667 

P2  =  13.33  psia 
CP/P0)1  =  0.1930 
(P/P0)2  =  0.7639 
PQ^  =  20.72  psia 
Po2  =  17.45  psia 
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